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We present the clustering properties of hard (2-8 keV) X-ray selected sources de- 
tected in a wide field (fa 2deg 2 ) shallow [fx (2 - 8keV) « 10 -14 ergcrn" 2 s _1 ] 
and contiguous XMM- Newton survey. We perform an angular correlation function 
analysis using a total of 171 sources to the above flux limit. We detect a ~ 4<j cor- 
relation signal out to 300arcsec with w(8 < 300 ) ~ 0.13 ±0.03. Modeling the two 
point correlation function as a power law of the form w(9) = (<9 O /0) 7_1 we find: 
d = 48.9^24 5 arcsec and 7 = 2.2 ± 0.30. Fixing the correlation function slope 
to 7 = 1.8 we obtain 80 = 22. 2^|'g arcsec. Using Limber's integral equation and 
a variety of possible luminosity functions of the hard X-ray population, we find a 
relatively large correlation length, ranging from r ~ 9 to 19 h~ 1 Mpc (for 7 = 1.8 
and the concordance cosmological model), with this range reflecting also different 
evolutionary models for the source luminosities and clustering characteristics. 



1. Introduction 

The overall knowledge of the AGN clustering using X-ray data comes mostly 
from the soft X-ray band (Boyle & Mo 1993; Vikhlinin & Forman 1995; 
Carrera et al. 1998; Akylas, Georgantopoulos, Plionis, 2000; Mullis 2002), 
which is however biased against absorbed AGNs. Recently, Yang et al. 
(2003) performing a counts-in cells analysis of a deep (f2-8keV ~3x 10~ 15 
erg s _1 cm -2 ) Chandra survey in the Lockman Hole North- West region, 
found that the hard band sources are highly clustered with ~ 60% of them 
being distributed in overdense regions. 

In this study we use a hard (2-8 keV) X-ray selected sample, compiled 
from a shallow (2-10 ksec per pointing) XMM- Yew ton survey near the North 
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and South Galactic Pole regions. A total of 18 pointings were observed out 
of which only 5 were discarded due to elevated particle background at the 
time of the observation. A full description of the data reduction, source 
detection and flux estimation are presented by Georgakakis et al. (2003). 
Here we just note that it comprises of 171 sources above the ha detection 
threshold to the limiting flux of ,/x(2 — 8kcV) « 10~ 14 ergs -1 cm~ 2 . 

2. Correlation function analysis 

In the present study we use as the estimator of the 2-point correlation 
function the following: w(9) — /(Ndd/Ndr) — 1, where Ndd and Ndr is 
the number of data-data and data-random pairs respectively at separations 
9 and 9 + d9. In the above relation / is the normalization factor / = 
2Nr/ {Nd — 1) where Nd and Nr are the total number of data and random 
points respectively. To account for the different source selection and edge 
effects, we have produced 100 Monte Carlo random realizations of the source 
distribution within the area of the survey by taking into account variations 
in sensitivity which might affect the correlation function estimate. Indeed, 
the flux threshold for detection depends on the off-axis angle from the center 
of each of the XMM-Newton pointings. Since the random catalogues must 
have the same selection effects as the real catalogue, sensitivity maps are 
used to discard random points in less sensitive areas (close to the edge of the 
pointings). This is accomplished, to the first approximation, by assigning a 
flux to the random points using the Baldi et al. (2002) 2-10 kcV log A— log S 
(after transforming to the 2-8 keV band assuming T = 1.7). If the flux of 
a random point is less than 5 times the local rms noise (assuming Poisson 
statistics for the background) the point is excluded from the random data- 
set. We note that the Baldi et al. (2002) log N — log S is in good agreement 
with the 2-8 keV number counts estimated in the present survey. 

The results of our analysis are shown in Figure 1, were the line corre- 
sponds to the best-fit power law model w(9) = (9 /9)' 1 ~ 1 using the standard 
X 2 minimization procedure in which each correlation point is weighted by 
its error. We find a statistically significant signal with w{9 < 300 ) ~ 
0.13 ± 0.03 at the 4.3a confidence level using Poissonian errors. The best 
fit clustering parameters are: 9 Q = 48.9^24:5 and 7 = 2.2 ± 0.30, where the 
errors correspond to la uncertainties. Fixing the correlation function slope 
to its nominal value, 7 = 1.8, we estimate 9 a = 22.2lg'g arcsec. 

Our results show that hard X-ray sources are strongly clustered, even 
more than the soft ones (see Vikhlinin & Forman 1995; Yang et al. 2003; 
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Figure 1. The 2-point angular correlation function of the hard (2-8 keV) X-ray sources. 
Insert: Iso-Ax 2 contours in the 7-S0 parameter space. 

Basilakos et al. in preparation). Our derived angular correlation length 
9 is in rough agreement, although somewhat smaller (within la) with the 
Chandra result of 9 = 40 ± 11 arcsec (Yang et al. 2003). The stronger 
angular clustering with respect to the soft sources could be either due to 
the higher flux limit of the hard XMM- Newton sample, resulting in the 
selection of relatively nearby sources, or could imply an association of our 
hard X-ray sources with high-density peaks 

3. The spatial correlation length using w{6) 

The angular correlation function w{9) can be obtained from the spatial 
one, £(r), through the Limber transformation (Peebles 1980). If the spatial 
correlation function is modeled as £(r, z) = (r/r )~ 7 (l + z)^^ +f \ The 
angular amplitude 9 is related to the correlation length r Q in three di- 
mensions via Limber's equation. Note that if e = 7 — 3, the clustering is 
constant in comoving coordinates (comoving clustering) while if e = —3 the 
clustering is constant in physical coordinates. We perform the Limber's 
inversion in the framework of the concordance ACDM cosmological model 
(On = 1 - fi A = 0.3, H = 70km s" 1 Mpc" 1 ). 

The expected redshift distribution and the predicted total number, N, 
of the X-ray sources which enters in Limber's integral equation can be 
found using the hard band luminosity functions of Ueda et al. (2003). 
We also use different models for the evolution of the hard X-ray sources: a 
pure luminosity evolution (PLE) or the more realistic luminosity dependent 
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density evolution (LDDE; Ueda et al 2003). The LDDE model with respect 
to the PLE gives an expected redshift distribution shifted to larger redshifts, 
with a median redshift of z ~ 0.75. 

For the comoving clustering model (e = 7 — 3) and using the LDDE 
evolution model, we estimate the hard X-ray source correlation length to 
be: r = 19 ± 3 /i" 1 Mpc and r = 13.5 ± 3 /tT 1 Mpc for 7 = 1.8 and 
7 = 2.2 respectively. While if e = — 3 the corresponding values are: r = 
11.5 ± 2 h^ 1 Mpc and r a = 6 ± 1.5 h^ 1 Mpc, respectively. 

The estimated clustering lengths (for 7 = 1.8) are a factor of > 2 larger 
than the corresponding values of the 2QZ QSO's (Croom et al. 2001). 
However, the most luminous, and thus nearer, 2QZ sub-sample (18.25 < 
bj < 19.80) has a larger correlation length (~ 8.5 ± 1.7 h^ 1 Mpc) than the 
overall sample (Croom et al. 2002), in marginal agreement with our e = — 3 
clustering evolution results. 

The large spatial clustering length of our hard X-ray sources can be 
compared with that of Extremely Red Objects and luminous radio sources 
(Roche, Dunlop & Almaini 2003; Overzier et al. 2003; Rottgering et al. 
2003) which are found to be in the range r ~ 12 — 15 h^ 1 Mpc. 
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